Linear kinetic Sunyaev-Zel'dovich effect and void models for acceleration 
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There has been considerable recent interest in cosmological models in which the current apparent 
acceleration is due to a very large local underdensity, or void, instead of some form of dark energy. 
Here we examine a new proposal to constrain such models using the linear kinetic Sunyaev-Zel'dovich 
(kSZ) effect due to structure within the void. The simplified "Hubble bubble" models previously 
studied appeared to predict far more kSZ power than is actually observed, independently of the 
details of the initial conditions and evolution of perturbations in such models. We show that the 
constraining power of the kSZ effect is considerably weakened (though still impressive) under a fully 
relativistic treatment of the problem, and point out several theoretical ambiguities and observational 
shortcomings which further qualify the results. Nevertheless, we conclude that a very large class of 
void models is ruled out by the combination of kSZ and other methods. 

PACS numbers: 98.80.Es, 95.36.+X, 98.65.Dx 



I. INTRODUCTION 

Probably the most surprizing cosmological discovery 
in recent decades has been the unexpected faintness of 
distant Type la supernovae (SNe) (see [1] for a review of 
the discovery). This has been widely interpreted as evi- 
dence for the accelerated expansion of the Universe due 
to a cosmological constant (or similarly behaving "dark 
energy"), or due to a modification to general relativity 
on very large scales. 

The standard picture of an accelerating, homo- 
geneous and isotropic Fricdmann-Lemaitre-Robertson- 
Walkcr (FLRW) Universe is now supported by a wide 
variety of observations in addition to SNe (see [2, 3] for 
reviews of the evidence). Nevertheless, uncertainty about 
the fundamental origin of a small cosmological constant 
has led to considerable effort in studying alternative mod- 
els in which inhomogeneity plays a key role. The simplest 
of these, both conceptually and practically, is the idea 
that an observer in an approximately spherical under- 
density, or void, which extends to rcdshift z ~ 1, will 
observe an apparent acceleration due to the inhomogene- 
ity [4-6] . These models dispense with the need for a mys- 
terious dark energy component at the cost of a violation 
of the Copernican principle, since the observer must be 
very close to the centre of the void in order to satisfy cos- 
mic microwave background (CMB) dipole constraints [7- 
10] . Although void models (and others based on inhomo- 
geneity) potentially solve the coincidence problem with 
the standard AFLRW model 1 by linking the appearance 
of acceleration to the formation of nonlinear structures, 
they are not free of temporal tuning [11]. 

Considerable effort has gone into confronting such void 
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1 We refer to the standard cosmological model as AFLRW to em- 
phasize its geometry, whereas the more common term "A cold 
dark matter" (ACDM) emphasizes its matter content. 



models for acceleration with a variety of different data 
(see, e.g., the reviews [12, 13] in this focus section). SN 
data can always be satisfied by choosing the radial pro- 
file of a void to match the redshift-luminosity distance 
relation of AFLRW [14]; nevertheless, these data do put 
crucial constraints on the depth and width of a void. 
Another class of observations involves the properties of 
structures within the void, e.g. studies of the baryon 
acoustic oscillations (BAO) [15, 16]. While potentially 
very constraining, these approaches are hampered by the 
difficulty of predicting the evolution of perturbations on 
LTB backgrounds [17, 18], although in special cases, such 
as near the centre, the evolution can be followed analyti- 
cally and strong constraints made [19]. More fundamen- 
tally, these approaches are hindered by our ignorance of 
the initial conditions (ICs) for the perturbations within 
the void: since the origin of the void itself is unclear, we 
cannot be sure that the ICs which apply at the radius of 
last scattering (and are manifest in the CMB) also apply 
locally. 

The primary anisotropics in the CMB contain a wealth 
of information, and hence their ability to constrain void 
models has been carefully studied. As first pointed out 
in [16], as long as the primordial spectrum is close to 
scale invariant, observations of the CMB temperature 
anisotropy spectrum can only be satisfied in void models 
if the local expansion rate is so low as to rule the models 
out. These results were shown to persist when the void 
is embedded in a spatially curved background [19-21]. 

The CMB is also expected to contain a variety of sec- 
ondary anisotropics generated at late times. Galaxy clus- 
ters within the void will see very large CMB dipoles and 
hence are expected to produce substantial anisotropics 
via the kinetic Sunyaev-Zel'dovich (kSZ) effect [22], as 
first suggested in [23]. This has been used to put con- 
straints on void models in [24, 25] . Recently, a related ap- 
proach in which the kSZ anisotropies due to all structure 
within the void are considered was introduced by Zhang 
and Stebbins [26] (hereafter ZS10). The authors found 
that void models predicted far more kSZ power than 



2 



is actually observed at small angular scales, and hence 
claimed that all such models were ruled out. However, 
ZS10 employed a simplified non-relativistic void model 
and did not examine the dependence of the effect on a 
variety of parameters. 

In this work we aim to repeat the analysis of ZSfO, 
but in a fully consistent general relativistic framework. 
We employ the exact, spherically symmetric pressureless 
matter solution to Einstein's equations, known as the 
Lemaitre-Tolman-Bondi (LTB) spacetime [27-29]. We 
assume that the LTB model contains no decaying mode, 
since we can then maintain the standard picture of an 
early homogeneous Universe consistent with inflation. 
We also ignore the possibility of a large, early isocur- 
vature mode between radiation and matter. We assume, 
with ZSf 0, that the matter power spectrum today is that 
of the standard AFLRW model, hence bypassing all un- 
certainty related to the ICs and evolution of perturba- 
tions inside the void. However, we show that the question 
of which scales the spectrum is to be evaluated at is not 
trivial. We confirm that void models which satisfy the 
SN data overpredict kSZ power, though considerably less 
so than indicated in ZS10, and point out that there are 
numerous theoretical ambiguities and observational un- 
certainties that affect the reliablity of kSZ calculations in 
void models. 

We begin in section II with a calculation of the kSZ 
power in LTB models. We point out that there is 
considerable ambiguity in relating length scales in LTB 
and standard AFLRW models, but propose a resolution. 
Next, in section III, we illustrate kSZ spectra and the 
distribution of power in wave number and redshift space. 
Section IV presents void model constraints using SN and 
local Hubble rate data. Section V examines the de- 
pendence of the kSZ power on various parameters. We 
extend our results to a calculation of the Compton y- 
distortion in section VI. Our conclusions are presented 
in section VII. The Appendix presents a brief descrip- 
tion of our LTB models. Throughout this paper we set 
c = I and define the conventional dimensionlcss Hubble 
rate ho via H = I00/iokms _1 Mpc _1 . 



II. CALCULATING THE LINEAR KINETIC 
SUNYAEV-ZEL'DOVICH EFFECT 

A. Derivation of the kSZ power 

The kSZ effect [22] due to Thomson scattering of the 
CMB from free electrons is conventionally written as 



AT(n) 
T 



J v(n, z) ■ n [1 



<5(n, z)]dr, 



(1) 



where n is the line-of-sight direction, S(n, z) = 
<5p e (n, z)j ' p e (z) is the comoving free electron density per- 
turbation in direction n and at redshift z, and r is the 
optical depth along the line of sight. The quantity v(n, z) 
is the relative velocity between the free electrons and the 



CMB rest frame (the frame in which the CMB dipolc 
vanishes) at (n, z). A scatterer at (n, z) will observe a 
radial component of a CMB dipolc directly related to the 
radial component /3(n, z) = v(n, z) ■ n. 2 

In an FLRW background, the dipole /3(n, z) vanishes 
by isotropy. Therefore in realistic models with structure, 
/3(n, z) is a perturbative quantity, corresponding to pe- 
culiar velocities. In an inhomogeneous LTB background, 
on the other hand, the lack of isotropy away from the 
centre means that the dipole will not generally vanish. 
Therefore, in a model with structure on top of an LTB 
background, we can decompose the dipole according to 
/3(n, z) = (3(z) + S(3(n, z), with the first part due to the 
LTB background (a very-large-scale "bulk velocity" ) , and 
the second part due to the superimposed perturbative 
structure (peculiar velocities relative to the LTB back- 
ground). Therefore we can write the kSZ anisotropy as 



= \p(z)+5p(n,z)][l + 5(n,z)]dT (2) 



AT(n) 
T 

(we henceforth drop the bar over f3(z)). We initially as- 
sume that the free electron fluctuations match those of 
the total matter, so that Sp e /p e = Sp m /p mi but consider 
the validity of this approximation in section VB. The 
dipole and density perturbations are evaluated on the 
light cone, so, e.g., 



5(n,z) = 5(n,t(z),r(z)) 



(3) 



for time coordinate t and radial coordinate r. We will 
often use the shorthand notation 



or 



Q(r) = Q(t(z),r(z)) = Q(z) 



Q(t) = Q(t(z),r(z)) 



(4) 



(5) 



for quantities Q(t,r) evaluated on the observer's past 
light cone. 

As we mentioned above, in FLRW backgrounds the 
background dipole vanishes, j3(z) = 0. In this case, the 
linear term in (2), J 5/3dr, is negligible due to a geomet- 
rical cancellation [30]. Therefore the second order term 
J SfiSdr dominates, and the effect in that case is some- 
times known as the Ostriker-Vishniac effect [31]. In LTB 
void models, the term J (3dr contributes only a small 
monopole 3 and hence we will ignore it. We will assume 
that the term J 5f3dr vanishes (or at least that it is sub- 
dominant) in LTB models as it does in the FLRW case. 
Hence we are left with one relevant first order term, 



AT(n) 
T 



(3(z)S(n, z)dr. 



(0) 



2 In terms of the spherical harmonic coefficients ai m to be defined 
in (13), we have /3(n, z) = *J 3/(4-7r)aio(n, z), when the polar 
axis is aligned along n. 

3 However, this term should lead to higher multipole anisotropics 
in nonspherical voids. 
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The fact that (6) is linear in the fluctuation amplitude 
leads us to call it the linear kSZ effect, after ZS10, and 
will result in much larger kSZ power in void models than 
in FLRW models. It is important to stress that the lin- 
ear kSZ effect is fundamentally different from the non- 
linear effect in FLRW models. Because the linear kSZ 
effect does not contain dipole perturbations (i.e. peculiar 
velocities), it is independent of the details of the veloc- 
ity power spectrum. On the other hand, the velocity 
power is crucial in the calculation of kSZ in FLRW mod- 
els. The insensitivity to the velocity power removes a 
significant source of uncertainty in the LTB kSZ calcu- 
lation, although considerable ambiguities will remain, as 
we will see. 

In order to evaluate the kSZ integral (6), we will wish 
to perform a harmonic decomposition of the perturbation 
field 8(n,z). To be able to do this, the field should be 
defined on a flat (or at least constant-curvature) space- 
like hypersurface. However, flat slices are not natural 
in general LTB spacetimes: they will not generally be 
orthogonal to the comoving worldlincs, so hypcrsurfaccs 
of constant proper time will not be of constant curva- 
ture. Nevertheless, recalling that we are restricting the 
LTB models to growing mode profiles, we can bypass this 
difficulty by relating the perturbation field on the past 
light cone, 5(n, z), to that on some slice t\, early enough 
that the LTB spacetime is close to FLRW at that time. 
We can then perform a decomposition on the early slice. 
Therefore, we introduce a growth function D(t) to relate 
the density perturbations at t- 1} <5i(n, r) = <5(n, U,r), to 
those at sonic later time via 



8(n,t(z),r(z)) = D(t(z))Si{n,r(z)). 



(7) 



By virtue of this expression, the coordinate r must be 
comoving. For the case of linear fluctuations about an 
FLRW background dominated by dust and cosmological 
constant at late times, we have the simple relation 



£>(*) = 



a(t) g(t) 

o(*i) g(U)'' 



(8) 



where a(t) is the FLRW scale factor and g(t) the usual 
growth suppression factor. 

It is very important to point out that the relation (7) 
neglects a couple of significant physical effects. First, 
in relating the perturbation at some coordinates (n, t, r) 
to that on the same comoving worldline at (n, t\,r), it 
ignores any scale dependence in the evolution. In par- 
ticular, at late times the nonlinear growth is expected to 
be greatest on the smallest scales. The second effect is 
that the scalar perturbation evolution on the LTB back- 
ground is expected to couple to vector and tensor modes, 
necessarily coupling different comoving worldlincs. How- 
ever, in the approximation that the coupling to tensors is 
ignored, it can be shown that the evolution satisfies a re- 
lation similar to (7) [17]. More importantly, as mentioned 
in the Introduction, we will not attempt a description of 
the perturbation ICs and evolution in this work; instead, 
we will assume that the actual matter power spectrum 
on the light cone is close to that of the standard AFLRW 
model. Therefore, we will consider it reasonable to ignore 
this latter effect. On the other hand, it will be necessary 
to deal with the former effect in an ad hoc manner by 
introducing scale dependence below. 

The LTB growing mode assumption allows us to ex- 
pand the field <5i(n,r) in spherical harmonic functions, 
>£ m (n), and spherical Bessel functions of the first kind, 
je(kr), according to 

<Si(n,r) = J- ( dk k^2S i: e m (k)ji{kr)Y em (n), (9) 

im 

where k is the wave number. It is important to stress 
the meaning of the quantities r and k here: the radial 
comoving coordinate r must be proportional to proper 
distance at ti (at least at background level, on the nearly 
FRLW background), and the wave number k is similarly 
proportional to proper wave number, in order that the 
harmonic decomposition be valid. 



Now, substituting (9) and (7), we can write (6) as 



AT(n) 



J- [ drF{r)D(t{r)) f dkkJ^Si 

^ Ira 



(k)jt(kr)Yt m (n), 



where 



dr 

F(r) = f3(r) — . 



(10) 



(11) 



As we mentioned above, this expression for the kSZ anisotropy ignores scale dependence in the perturbation evolution. 
However, as we will see below, the matter power on nonlinear scales will be crucial to the strength of the kSZ technique 
in constraining LTB models. Therefore, in order to capture the effect of the nonlinear power, we must promote the 
growth function to be scale dependent, D(t) — > D(t, k). With this replacement, (10) becomes 



= \/f / drF (r) J dkkD(t(r),k)J2Si 



( in 

(k)j e (kr)Y em (n). 



(12) 
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To avoid this ad hoc procedure would require developing 
the theory of perturbations on LTB backgrounds. While 
some progress has been made in the linear case [17, 18], 
the nonlinear case has not yet been addressed. However, 
although strictly the growth function D(t, k) is ill-defined 
since we cannot perform a harmonic decomposition at 
late times, this procedure might be justified in that it 
may be possible to consider such a decomposition valid 
on the very small scales (much smaller than the LTB 
curvature scale) relevant to the kSZ effect. This leads to 
our first major caveat: our result will likely be invalid 
on the largest angular scales, and the degree of accuracy 
will be uncertain even on small scales. 

Next we can calculate the spherical harmonic coeffi- 
cients of the observed temperature field, 



_ f AT(n) 



T 



Y; m (n)dQ 



(13) 



drF(r) / dkkD{r,k)5- htm {k)u{kr). (14) 



(15) 



Finally, the kSZ power at multipole i is given by 
Ci = (a em a* im ) 

= 4tt J yPftCfc) (y drF(r)D(r,k)h(kr)j . (16) 

In deriving this final expression, we have used the defin- 
ing relation 

(8 Um {k)5l P , ml (k')) = ^-V Si (k)5(k - k')8u,8 mm , (17) 

for the power spectrum T > s i (k) of Gaussian random mat- 
ter field <5i(n, r). The spectrum Vs^k) is well defined 
since a harmonic decomposition is possible at t\, when 
the spacetime is near FLRW. 

Equation (16) provides a general expression for the 
kSZ power on LTB backgrounds, subject to the major 
caveat regarding the ambiguity of defining the growth 
factor D(r,k). The factor F(r) in the integrand will de- 
pend on the particular radial LTB profile chosen. It can 
be written 



where 



and 



F(r) = /?(r)^, 
dz dr 



dr_ _ OT/b(2 - Yn c )p m (z) 
dz ~ 2m p (l + z)H\\{z) 



dz , , , s Y' 



(18) 



(19) 



(20) 



Here or is the Thomson cross section, /b = pb/Pm is 
the baryon fraction, Yhc is the helium mass fraction, m p 
is the proton mass, and the remaining LTB functions 



are defined in the Appendix. This expression applies af- 
ter reionization, when we have assumed that the bary- 
onic matter is completely ionized. We use the value 
Yhc = 0.24 throughout this work, and nominally set 
fb = 0.168, as implied by CMB observations [32]. As 
we discuss in detail in section VA, it is likely that this 
standard value of baryon fraction should be modified lo- 
cally in void models, leading to considerable uncertainty 
in our kSZ predictions. 



B. The matter power spectrum 

After the LTB profile is specified, the remaining quan- 
tities required to calculate the kSZ power are the matter 
power spectrum Vs x {k) and the growth function D(r, k). 
As explained above, we do not attempt to calculate 
these quantities from ICs; rather, we begin by assuming 
that the actual matter power along our past light cone 
matches that of the standard AFLRW model. Later we 
will consider the effect of relaxing this assumption. How- 
ever, implementing even this simple prescription for the 
matter power is far from trivial. First, as already men- 
tioned in section II A, it is not clear how to define a power 
spectrum at late times when a harmonic decomposition is 
not possible. But there is a second important reason: the 
wave numbers k with respect to which the power spectra 
are specified have different meanings and behaviours in 
LTB and FLRW models. Since the matter power spec- 
trum can depend sensitively on scale, this will lead to 
a further significant ambiguity in specifying the matter 
power. 

To see this, note first of all that, for FLRW models, 
the power spectrum on the past light cone, V${k,z), is 
usually specified in terms of comoving wave number k (or 
k/ho), which is equivalent to proper wave number today. 
In general, perturbation modes in LTB and FLRW mod- 
els with identical proper wave numbers today will not 
share identical proper wave numbers on the past light 
cone, due to the very different background evolution of 
the two models. (Recall that the LTB and standard 
AFLRW background evolutions differ at zeroth order.) 
Therefore it is certainly incorrect to state that the LTB 
model should share the same Vs{k, z) at the same proper 
wave number today as the standard AFLRW model. 

However, it is also generally incorrect to require that 
the power Vs(k,z) be the same in both models at the 
same proper k specified at the same redshift z on the past 
light cone. The reason is that what is actually observed 
in galaxy surveys is not proper wave numbers at some 
z, but rather angular and redshift separations. A pair of 
galaxies with the same proper separation at the same z in 
LTB and FLRW models will generally be observed with 
very different angular and redshift separations in the two 
models, again due to the different background geometry. 

Therefore, what we require is that the LTB model has 
the same matter power on the same angular and redshift 
scales, at the same z, as the standard AFLRW model. In 
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any model, for a mode oriented perpendicular to the line 
of sight with proper wave number k± at redshift z, the 
corresponding angular scale on the sky is 



AO = 



k±(z)d,A(z) 



(21) 



where oIa(z) is the angular diameter distance to z, mea- 
sured as a proper distance at z, and a is a constant (in- 
dependent of the model) of order unity. Therefore, if the 
two models are to share the same power at the same an- 
gular scales, they must have the same power at proper 
wave numbers related by 



^LTB / \ _ j A / \ d A( z ) 
"± \ z ) ~ K ±\ z ) iLTB 



(22) 



In practice, if the LTB model is to fit the Type la super- 
nova data and measurements of the local Hubble rate, 
then d A TB (z) must be similar to d A (z), at least out to 
z ~ 1, and hence a first approximation would be to set 
k LTB {z) = k A {z). 

Similarly, if the two models are to share the same power 
at the same redshift scales, they must have the same 
power at proper wave numbers related by 



LTB 



fcf(z)- 



LTB 



(z) 



H A (z) 



(23) 



This leads directly to a fundamental problem: generi- 
cally, we will have fc LTB (z) ^ fc LTB (z), when k A (z) = 

k A (z), due to the background geometry of the LTB mod- 
els. 4 This means that there is no unique scale in AFLRW 
corresponding to a particular scale in the LTB model, 
and vice versa. (Indeed, this discrepancy between angu- 
lar and redshift scales in the two models results in the 
strength of the radial BAO scale in constraining LTB 
models [16, 19].) 

To patch over this problem, we define an "isotropized" 
version of the relation between wave numbers in the two 
models by 



k A (z)^ (k A (z)Yk A (z) 



1/3 



,LTB 



d A (z) 



■LTB 



CO 



H A (z) 
Hf B (z) 



(24) 

1/3 

(25) 



This relation weights the transverse relation more heav- 
ily than the radial one, since there are twice as many 
transverse dimensions. It is important to stress that the 
particular form of this relation we have chosen is still 
largely arbitrary. [This notion of an isotropized scale is 



4 This same model dependence in the meaning of Vs(k, z) will also 
be present when comparing different FLRW models, although to 
a lesser degree than with LTB models. 



akin to the isotropized distance measure often used in 
the literature on BAO observations (e.g. [33])]. Our fi- 
nal statement, then, is that the LTB and AFLRW mod- 
els must share the same matter power on scales related 
by (25). The ambiguity in this relation leads to our sec- 
ond major caveat: uncertainty in the k scales, together 
with the sensitivity of the matter power to k in standard 
spectra, will lead to further uncertainty in the calculated 
kSZ power. 

The evolution of proper wave numbers in the radial 
and transverse directions in the LTB model is straight- 
forward to specify. For a mode of proper wave number 
k(U) specified at the time t\ early enough that the space- 
time is near FLRW (and hence the power spectrum is 
presumably isotropic), we have 



kj_(t,r) = k(ti,r) 



Y(t u r) 
Y(t,r) ' 



k\\(t,r) = k(t u r) 



Y'(t u r) 



Y'(t,r) 
(26) 

We stress again that our inability to perform proper har- 
monic decompositions of fields at late times implies that 
we must interpret relations of this sort loosely. 

Recall that the general expression for kSZ power (16) 
involves an integral over proper k at t;. To determine the 
matter power on the corresponding scale, we first use (26) 
to relate the k values at ti to those on the light cone in 
the LTB model. We then use (25) to calculate the cor- 
responding scale in the standard AFLRW model, using 
the fitting function of [34] and the WMAP7 best-fit pa- 
rameters [32] to evaluate d^(z) and H A (z). The AFLRW 
power can be calculated using any of several public codes; 
we used CAMB [35], which optionally includes nonlincar 
power with a Halofit calculation [36]. Our matter spec- 
trum used the normalization provided by CMB observa- 
tions [32], us = 0.81, where as is the linear amplitude 
on a scale of 8/ ho Mpc. As we discuss in section VA, 
however, there is considerable ambiguity over this choice 
in the context of LTB models. Finally, note that (25) 
becomes ill-defined when _ff LTB (z) < 0. This does in 
fact happen for the very deepest voids, leading to multi- 
valued distance- redshift relations [37]. Therefore our kSZ 
calculations will not be valid in the multi- valued regime. 



C. The Limber approximation 

We can significantly simplify the kSZ power expres- 
sion (16) using the so-called Limber approximation [38]. 
This approximation is purely mathematical, as opposed 
to the physical approximations discussed previously To 
employ it, we substitute 



je(kr) 



2i + l 



+ 1/2- kr) 



(27) 



in (16). This substitution will be accurate to (D(l/£ 2 ) 
when the integrand F(r)D(r, k) is slowly varying on the 
scale of the Bessel oscillations [39]. The result for the 
kSZ power is 



6 



Of 



4tt 2 



2£+l 
16tt 2 



dk 



21+1 
2k 



D 



(2£+l) 



— J drrF 2 {r)D 2 



21 



21+1 
2k 
1 



2r 



k 2 

'21+1 
2r 



(28) 
(29) 



If we attempt to define a matter power spectrum on the 
light cone by T s {k,z) = D 2 (z,k)T8 1 (k), then (29) be- 
comes 

c <=(^/*^< r K^ w )- <30) 

Note that all k values in these expressions arc proper 
wave numbers evaluated at To specify the standard 
AFLRW matter spectrum, we must relate these early 
time k values to corresponding scales in the AFLRW 
model using the procedure described in section II B. 

Equation (30) apparently has a simple interpretation 
as an integration down the past light cone of a weighted 
local matter power spectrum. However, it is important 
to stress that the interpretation of the late time spec- 
trum Vsik^z) is highly ambiguous, as discussed in sec- 
tions II A and II B. The same replacement of Vs^k) with 
Vs(k,z) cannot be made with the exact expression for 
kSZ power, (16), and hence neither can such a simple 
interpretation be made. Note that (30) is equivalent to 
the corresponding expression in ZS10, up to 0(l/£). The 
specific ^-dependent prefactors in (30) ensure that it is 
accurate to 0{l/£ 2 ). 

The Limber approximation provides a significant ad- 
vantage to numerical computations, especially if a large 
region of parameter space is to be explored. Therefore 
it is important to check its accuracy. Figure 1 presents 
the relative error between the Limber calculation of kSZ 
power using (29) and the exact calculation using (16), 

A /^Limber v^cxact 

/^•Gxact /^cxact ' ^ ' 

for the fiducial LTB model (which is described below). 
The Limber approximation can be seen to be very ac- 
curate, to better than a tenth of a percent, over a wide 
range of angular scales. This was expected, as the inte- 
grand F(r)D(r, k) in (16) varies slowly over the width of a 
void. 5 Henceforth, all calculations will be performed us- 
ing the Limber approximation, (29), or equivalently (30). 



° A more rapidly varying LTB profile might result in a less accurate 
Limber approximation. 
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FIG. 1: The relative error of the Limber approximation given 
by (29), C^ lmber , with respect to the exact kSZ power given 
by (16), C*r act . ACf = Ct imbcr - CT act . Results are pre- 
sented for the fiducial LTB model. 



III. BEHAVIOUR OF KSZ POWER IN VOID 
MODELS 

With the formalism for calculating the kSZ power in 
place, we can now procede to calculate spectra for spe- 
cific LTB models, and to investigate where in k and z 
space the kSZ power is sourced. First, the LTB back- 
ground must be specified. This entails the specification 
of a radial profile K(r), together with the local Hubble 
rate at the centre today, Hq. Our fiducial profile for K(r) 
is defined in the Appendix. Our fiducial model is further 
specified by the parameters fij°^ = 0.2, zl = 0.5, and 
h = 0.71, where n } ^ = il\° c (z = 0) is the local den- 
sity parameter at the centre today, and zl is a measure 
of the width of the profile in redshift. As we will see in 
section IV, the fiducial model is a very good fit to the 
supernova data. 

Then the kSZ power can be calculated using the Lim- 
ber expressions (29) or (30). Importantly, in evaluat- 
ing the function F(r) using expressions (18) to (20), all 
background quantities must be evaluated consistently us- 
ing the specified LTB background. The dipole (3{r{z)) is 
also calculated consistently in the LTB model by propa- 
gating past-directed null rays radially inwards and out- 
wards from the point (t(z),r(z)) on the past light cone 
of the central observer to an early time, and comparing 
the resulting redshifts as described in [19]. The matter 
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FIG. 2: The kSZ power calculated in the Limber approxima- 
tion using the Halofit nonlinear matter power (dashed, red 
curve) and linear power (dotted, blue curve), compared with 
the temperature anisotropy power in the standard AFLRW 
model (solid, black curve). Results are presented for the fidu- 
cial LTB model. 
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FIG. 3: The derivative of the kSZ power with respect to 
redshift, calculated in the Limber approximation using the 
Halofit nonlinear power for three multipoles: I — 30 (dot- 
ted, blue curve), I = 300 (dashed, red curve), and i — 3000 
(solid, black curve). Results are presented for the fiducial 
LTB model. 



power spectrum in the LTB model is assumed to match 
that of the standard AFLRW model, and is calculated as 
described in section II B. 

A comparison of the kSZ power spectrum calculated for 
the fiducial LTB model with the temperature anisotropy 
power in the standard AFLRW model is presented in 
figure 2. The quantities plotted are the conventionally 
scaled power spectra, 



2ir 



(32) 



Figure 2 shows the kSZ power calculated using both the 
nonlinear (Halofit) matter power as well as the corre- 
sponding linear spectrum. Nonlinearities in perturbation 
evolution enhance the matter power on small scales. For 
the fiducial LTB model and the nonlinear matter power, 
the kSZ power is found to dominate over the standard 
temperature anisotropics for I > 1000. Clearly the pres- 
ence of such a large kSZ component would result in a 
total anisotropy spectrum drastically different from that 
actually observed. Although the kSZ power is reduced 
considerably in the linear matter power case, it still ex- 
ceeds the standard temperature anisotropies by a large 
factor for the largest multipoles. While the fiducial model 
would therefore appear to be immediately ruled out, we 
will see in section VA that the kSZ power depends sen- 
sitively on several parameters. 

Figure 2 suggests that the nonlinear matter power on 
small scales is the main source of the apparently very 
strong constraint that the kSZ can provide on void mod- 
els. We can examine this more explicitly by plotting 
derivatives of the kSZ power. Figure 3 shows the deriva- 
tive with respect to redshift, 



dV e _ S^(£ + 1) dr 2 
dz ~ (2£+l)3 dz 1 >ld 



2£+l 
2r ' 



z(r) . (33) 



In words, it shows how the contributions to the total 
power at various distances along the line of sight are 
distributed. The figure shows that, essentially indepen- 
dently of £, most of the kSZ power originates at redshifts 
z ~ 0.5, i.e. near the edge of the void for the fiducial 
model which extends to Zi = 0.5. The redshift distribu- 
tion of the kSZ power is largely determined by the dipole 
function, (3(z). Indeed, the dipole peaks near z = zl and 
exhibits a zero at z ~ 2 which is visible in figure 3. 

Figure 4 shows the derivative of kSZ power with re- 
spect to logarithmic k interval, 



dV t 2tt£(£ + 1) 2 
- F 



dhi k 



21+1 



2£+l\ T s (k,z) 



2k 



k 2 



(34) 



for both the nonlinear and linear matter power spectrum. 
That is, it shows how the contributions to the kSZ power 
from different length scales are distributed. Importantly, 
the distributions are plotted with respect to the scale k A , 
calculated as described in section II B (and measured as 
a proper wave number today). This is the scale that we 
would conclude a feature lies at if we mistakenly assumed 
that the AFLRW model was correct, when actually the 
specified LTB model was correct. 

As expected from the relation (21), figure 4 shows that 
the wave numbers sourcing the kSZ power scale in inverse 
proportion to the angular scale examined on the sky. In 
addition, it is clear that the length scales sourcing the 
most powerful part of the kSZ spectrum, near I ~ 3000, 
are strongly nonlinear. Importantly, the figure also shows 
that the length scales we would infer to make the largest 
contribution to the high-£ kSZ effect are very small: for 
I = 3000, the corresponding length scales are roughly 
1/3 Mpc. We will discuss the significance of these small 
scales in section VB. 
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FIG. 4: The derivative of the kSZ power with respect to wave 
number, calculated in the Limber approximation using the 
Halofit nonlinear power (top, solid curve of each pair) and lin- 
ear power (lower, dotted curves). Results are shown for three 
multipoles: I = 30 (leftmost, blue pair of curves), £ — 300 
(central, red curves), and £ — 3000 (rightmost, black curves). 
The distributions are plotted with respect to the scale k A /ho, 
calculated as described in section II B, and measured as a 
proper wave number today. Results are presented for the 
fiducial LTB model. 



IV. CONSTRAINTS ON VOID MODELS 

The results of the previous section show that LTB 
models can potentially exhibit very high kSZ power. In 
order to determine whether viable void models predict 
too much kSZ power, we must select models which fit 
other observations. For our analysis we used the minimal 
amount of data required to fix the void profile, and hence 
obtained the most independent kSZ constraints possible. 
To do this we used redshift-magnitude measurements of 
Type la SNe from the Union2 compilation [40], consist- 
ing of 557 SNe in the range z = 0.015-1.4. In our fitting 
we adopted the standard procedure of marginalizing over 
the unknown absolute magnitude, which is equivalent to 
marginalizing over Hq . A further input is therefore a local 
estimate of Hq (i.e. for z < 0.1), which is essentially free 
of any assumptions regarding the cosmological model. 
For this we used a prior of ho = 0.72±0.08 from the Hub- 
ble Space Telescope (HST) Key Project [41]. Since our 
choice of prior has a larger error bar than the more recent 
measurements in [42, 43], it is more conservative. Note 
that the SN data are not model independent; instead, 
they assume a standard AFLRW background (e.g., see 
the discussion in [44]). This will result in further uncer- 
tainty in the kSZ power, considering the strong sensitiv- 
ity of the power to the void width which we will explicate 
in section V A. 

In our previous analyses, we found that much lower 
values of the local Hubble rate (h ~ 0.45) were required 
for void models to fit the CMB primary anisotropics, as- 
suming a ncar-power-law primordial spectrum [16, 19]. 



Although such low values of Ho are in severe conflict with 
the local measurements in [41-43], this tension might in 
principle be circumvented by substantial fine tuning of 
the primordial spectrum, allowing for higher values of 
Ho ■ For this reason we chose a prior on Hq which is con- 
sistent with local observations, and we did not include 
the CMB primaries in our constraints. Another reason 
we excluded CMB primaries is that, as figure 2 shows, a 
large linear kSZ component will substantially modify the 
observed anisotropics over a wide range of scales. 

We used CosmoMC [45] to generate Markov-Chain- 
Monte-Carlo chains to estimate confidence limits. The 
LTB profile we used was the fiducial profile defined in 
the Appendix. Three parameters were varied: the void 
width L and depth Kq, along with Ho- For each sam- 
ple of the chain we computed several derived parameters: 
the redshift to r = L, zl, the central matter density pa- 
rameter, O^ , and the kSZ power, V e , at I = 2500, 3000, 
3500, 4000, and 4500. As in our previous analyses we ap- 
plied a conservative prior of O^o > 0.1 [46, 47]. Deeper 
voids, with fl^o < 0.1, must be correspondingly wider 
to fit the SN data and hence have the largest kSZ power 
(as we will see in figure 6). Hence the deepest voids are 
most strongly ruled out by the kSZ constraint. We also 
applied a prior of T 2 T)g < 10 000 /iK 2 , since such models 
would already be excluded at high significance by small- 
scale CMB experiments such as ACBAR [48]. Since we 
did not use primary CMB data, which fixes the baryon 
fraction /b (at least at the radius of last scattering) , we 
needed to choose a (local) value; we set the fiducial value 
of /b = 0.168 for the kSZ calculation. Similarly, we set 
the fiducial value of a a = 0.81 for the matter power am- 
plitude today. The validity of these assumptions is dis- 
cussed in the next section. 

The best-fit model has a total \ 2 = 539, and therefore 
a goodness-of-fit comparable to AFLRW and the void 
models in our previous work. In figure 5 we show a his- 
togram of the kSZ power from the chains at various £ 
values. Each distribution is similar over the £ range we 
consider, and hence lower limits on the kSZ power are 
largely independent of £. The maximum-likelihood value 
is approximately 2500 /iK 2 , with an extended tail due 
to deeper, wider void profiles. Asymmetric distributions 
make estimating confidence limits (from counting sam- 
ples in the histogram) dependent on the tail, but we have 
checked the robustness of our constraints by also exclud- 
ing samples with T^Vg > 5000 /iK 2 . The 2cr lower limit 
on the kSZ power at £ = 3000 is reduced from 1630 /iK 2 
to 1440 /xK 2 when applying this constraint. We take the 
lower (and more conservative) limit as our baseline value 
to compare with recent CMB observations on very small 
scales from the South Pole Telescope (SPT) [49, 50] and 
the Atacama Cosmology Telescope (ACT) [51]. In fig- 
ure 6 we show the 1,2, and 3cr confidence levels in the Zl- 
tt\° c plane. Note that our fiducial model, with zl = 0.5 
and ri|n C o = 0.2, is very close to the best fit model. 

Single frequency bandpowers from the 150 GHz and 
and 148 GHz channels of SPT and ACT limit the total 



9 




0.0 2000.0 4000.0 6000.0 8000.0 10000.0 
^[/xK 2 ] 

FIG. 5: Histograms of the kSZ power, calculated in the Lim- 
ber approximation using the Halofit nonlinear power for five 
multipoles: t = 2500 (solid, black curve), £ = 3000 (long 
dashed, red curve), £ = 3500 (short dashed, green curve), 
£ = 4000 (dotted, blue curve), and £ = 4500 (dot-dashed, or- 
ange curve). Results are presented for the fiducial LTB pro- 
file, fitted to the Union2 SN data, and using the HST prior 
of ho — 0.72 ± 0.08. Each curve is normalized to unity at its 
maximum. 



power (including primary CMB, point sources, thermal, 
and kinetic SZ) at I = 3000 to < 50 ^K 2 . While both 
SPT and ACT also present considerably tighter upper 
limits on the kinetic SZ component itself (and these lim- 
its were used in ZS10), such limits are model dependent. 
This is because the contribution of primary CMB in par- 
ticular may be quite different in LTB models than in 
AFLRW, recalling that the CMB and Hq cannot be si- 
multaneously fit without substantially modifying the pri- 
mordial spectrum. Therefore we consider the SPT and 
ACT total power to be a conservative upper limit to the 
kSZ power from voids. 6 Even taking this conservative 
approach, the predicted kSZ signal is at least a factor of 
roughly 30 higher than the observational upper limits, 
and hence rules out void models which fit the SN data 
at high significance when the fiducial parameters are as- 
sumed. 



V. ROBUSTNESS OF KSZ CONSTRAINTS 

A. Parameter dependencies 

Although the results in section IV suggest that the kSZ 
anisotropy power is large enough to immediately rule out 



The total powers will almost certainly be overestimates of the 
kSZ component, since they include a frequency-dependent com- 
ponent presumably due to point sources. 
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FIG. 6: Solid contours are the kSZ power, T Q 2 T>zooo (in fjK 2 ), 
calculated in the Limber approximation using the Halofit mat- 
ter power, as a function of central matter density parameter, 
fi™ , and a measure of the width of the void, zl- Results are 
presented for the fiducial LTB profile with ho = 0.71. Dotted 
contours are the 1, 2, and 3a confidence levels from the SN 
data, which used the HST prior of h = 0.72 ± 0.08. 



void models for acceleration, it is important to examine 
the robustness of the kSZ power as a test of homogeneity. 
In particular, how sensitive are the predictions to vari- 
ous model parameters? This was already addressed to 
some extent implicitly in section IV, where the LTB pro- 
file (width and depth), as well as local Hubble rate, were 
varied. In this section we will illustrate these dependen- 
cies explicitly, and examine the effects of the remaining 
parameters. 

Figure 6 shows the kSZ power for I — 3000 as a func- 
tion of a measure of the width of the void, zl, and the 
depth of the void as measured by the central density pa- 
rameter, fij^ c . The LTB profile used was the fiducial 
profile, with local Hubble rate ho = 0.71. Superimposed 
on figure 6 are the 1,2, and 3er confidence levels from the 
Union2 SN data calculated in section IV. The kSZ power 
is seen to be a very sensitive function of the void width, 
varying by roughly a decade over the region allowed by 
the SNe. Note that, for the fiducial profile, the multi- 
valued region of parameter space (for which H\\{z) < 0) 
occurs only for the deepest voids, with ^l\°-^ Q < 0.05. 
Therefore, as noted in section II B, the kSZ calculation 
is not valid in this regime and we do not plot it in fig- 
ure 6. (The excluded shell-crossing region occurs for even 
deeper voids than the multi- valued models.) 

We can understand the general forms of the parame- 
ter dependencies illustrated in figure 6 quite easily. For 
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fixed void width zl (and fixed Hq), the kSZ power in- 
creases as fi^ c decreases from unity. This is mainly due 
to the dipolc factor /3 2 in the expression for the kSZ, (29). 
Inside the void, we can define an effective curvature pa- 
rameter by = 1 — f^ c - Then, for a growing mode 
background solution, we expect Off ^ SH/H, where 
5H is some measure of the perturbation in expansion 
rate inside the void relative to the outside [17]. Finally, 
we expect the dipole at fixed rcdshift to be proportional 
to the expansion perturbation, 

/3 ~ DSH, (35) 

where D is some proper distance measure [52]. Combin- 
ing these relations we find 

Q cx /? 2 (z L ) cx (l-fi£?o) 2 - ( 36 ) 
Although there are considerable ambiguities in defining 
most of the quantities used to derive this relationship, 
it is satisfied reasonably well by the numerical results of 
figure 6. 

Next, for fixed ti^o (and fixed Ho), figure 6 shows 
that the kSZ power grows rapidly with increasing void 
width zl. Using the relation (35), we have (3 2 {r) cx r 2 , 
since coordinate r will be approximately proportional to 
a proper distance measure. Then, ignoring the scale de- 
pendence of the matter power spectrum, the kSZ power 
expression (29) implies 

C t oc [ r 3 dr cx L 4 cx z%, (37) 
Jo 

for sufficiently small zl- Again, this relation is indeed 
roughly satisfied by the numerical results. 

Figure 7 illustrates the effect of varying the local Hub- 
ble rate, H , while keeping the LTB depth (as measured 
by Q[£ c o) fixed. Larger H leads to larger kSZ power, and 
this dependence is due to a combination of background 
factors in the relations (18) to (20), together with the 
scale dependence of the matter power. The variation in 
power over our chosen prior range of ho = 0.72 ±0.08 [41] 
is not large, and of course the variation is consider- 
ably smaller over the range of the newer measurement, 
h = 0.738 ± 0.024 [43]. 

The remaining parameters that affect the kSZ 
anisotropy power are the linear matter power amplitude 
today, erg, conventionally expressed on a scale of 8 /ho 
Mpc, and the baryon fraction fy,. The fluctuation am- 
plitude could also be expressed in terms of the primor- 
dial perturbation amplitude, although in our approach 
of specifying the matter power at late times rather than 
evolving perturbation ICs, the parameter as is more nat- 
ural. If the fluctuations were linear, the kSZ power would 
simply scale as erf. However, in the nonlinear regime rel- 
evant to the kSZ effect for high £, we expect [53] 7 

C e oc Vs oc crg~ 3 . (38) 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



FIG. 7: The kSZ power, r 2 30 oo (in ^K 2 ), calculated in the 
Limber approximation using the Halofit matter power, as a 
function of local Hubble rate, -ffo, and a measure of the width 
of the void, zl- Results are presented for the fiducial LTB 
profile with Q}^ Q = 0.2. 



This means that, e.g., an uncertainty of 30% in erg would 
lead to an uncertainty in the kSZ power of roughly a 
factor of two. 

Of course, observations of the CMB are usually con- 
sidered to constrain er 8 very tightly, to a few percent [32] . 
But this is based on the assumption of homogeneity: void 
models, on the other hand, contain an extremely large 
inhomogencity of mysterious origin, and hence there is 
no necessary link between the perturbation amplitude 
locally and that at the last scattering radius. In addi- 
tion, as mentioned in the Introduction, if we are to sat- 
isfy local measurements of H n , the primordial spectrum 
must be modified substantially [16, 19], further clouding 
the connection between the local and CMB amplitudes. 
Local measurements of erg (or, more precisely, measure- 
ments at z ~ 0.5) are what we need to test void models 
via the kSZ effect; unfortunately, they are considerably 
more uncertain than those based on the CMB. For exam- 
ple, measurements based on weak lensing (e.g., [54]) con- 
strain only a combination of a% and Q m . Extracting con- 
straints on erg from these results would likely be strongly 
model dependent. Recalling that our calculations have 
used the standard CMB-based value, erg — 0.81 [32], we 



7 For the ordinary (nonlinear) kSZ effect in standard AFLRW mod- 



els, the kSZ power scales like an even larger power of erg, since 
the same amplitude also determines the velocity perturbation 
spectrum [53]. 
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can expect the uncertain local matter normalization to 
add significant ambiguity to our calculations. 

The final parameter determining the kSZ power is the 
baryon fraction /b- This affects the kSZ via the back- 
ground optical depth, which is related to the background 
electron density (assumed equal to the baryon density) 
through (19). A constant baryon fraction affects the kSZ 
power through the simple scaling 

Ct oc fl. (39) 

As was the case with the matter power amplitude, con- 
straints on /b from the CMB will not generally be ap- 
plicable in an inhomogeneous Universe, and, again, local 
measurements are not very constraining. In addition, in 
void models, the presence of cosmological-scale adiabatic 
(curvature) inhomogencity of mysterious origin suggests 
that we should also consider the possibility of isocurva- 
ture inhomogeneity between baryons and dark matter, in 
the form of an r-dependent /b [55]. Thus, recalling that 
in our calculations we have used the CMB-based value 
/b = 0.168 [32], our poor knowledge of the local baryon 
fraction in the context of LTB cosmologies will result in 
considerable further uncertainty in the kSZ power calcu- 
lations presented here. 

B. Small-scale baryonic astrophysics 

Recall from figure 4 that, for the large £'s which provide 
the strongest constraints on void models, the dominant 
length scales sourcing the kSZ in the fiducial model are 
very small, namely less than a Mpc. Matter perturba- 
tions on such small scales are expected to be strongly 
nonlinear, as we have already seen in section III when 
comparing the kSZ power calculated for linear and non- 
linear matter power spectra. However, when we began 
the derivation of the kSZ effect, we assumed that the 
free electron fluctuations match those of the total mat- 
ter, so that Sp e /p e = 5p m / 'p m . In the nonlinear regime, 
the baryonic (and hence electron) fluctuations are ex- 
pected to be suppressed with respect to the total (or 
dark) matter power, due to interactions that become im- 
portant on small scales (e.g., see [56]). Therefore, we 
expect the kSZ power sourced at these scales to be some- 
what smaller than the calculations thus far (which used 
the standard AFLRW total matter power spectra) have 
indicated. To help examine the importance of this ef- 
fect, it will be useful to consider the dependence of the 
sourcing length scales on the various model parameters. 

In figure 8 we illustrate the dependence on the width 
and depth of the void of the wave number providing the 
peak contribution to the £ = 3000 kSZ power, again for 
the fiducial profile and h = 0.71. (The plotted peak 
wave number maximizes dT>g/ 'dink, given in (34).) As 
we did for figure 4, we here plot the k scale that we would 
conclude sources the kSZ, if we mistakenly assumed that 
the correct model was standard AFLRW, when actu- 
ally the (fiducial) void model was correct. This effective 




FIG. 8: Solid contours are the peak wave number contribut- 
ing to the kSZ power at t = 3000, calculated in the Limber 
approximation using the Halofit nonlinear power, as a func- 
tion of central matter density parameter, f2^ , and a mea- 
sure of the width of the void, zl- The contour values are the 
peak scale k A /ho (in Mpc -1 ), calculated as described in sec- 
tion II B, and measured as a proper wave number today. Re- 
sults are presented for the fiducial LTB profile with ho = 0.71. 
Dotted contours are the 1, 2, and 3er confidence levels from 
the SN data, which used the HST prior of h = 0.72 ± 0.08. 



AFLRW scale is calculated according to the prescription 
of section II B. Figure 8 shows that the peak k scale gen- 
erally increases as the void width decreases. This can 
be understood simply from the geometrical relation (21): 
keeping £ (or the angular scale) fixed, the relevant k scale 
must be inversely proportional to z, since dx{z) oc z, 
for sufficiently small rcdshift. This simple relationship is 
modified for large z or very deep voids, since rclativistic 
effects then become important. 

For the entire region of parameter space illustrated in 
figure 8, the peak scale is smaller than 1 Mpc" 1 , and for 
the narrowest voids consistent with the SN constraints 
from section IV, namely zl — 0.35 at 2cr, the peak scale 
is k A /h ~ 4Mpc _1 . On such small scales, hydrody- 
namical simulations indicate considerable suppression of 
baryon power with respect to dark matter power, al- 
though the various studies differ quantitatively in their 
predictions (e.g., see [56, 57]). Observations, on the other 
hand, most directly probe either the galaxy or the total 
matter power, rather than the baryon power. Therefore, 
it is clear that the poorly understood details of small- 
scale baryonic astrophysics add further considerable un- 
certainty to our estimates of kSZ power, although in this 
case we can say that our calculations will overestimate 
the kSZ for large £'s. It is also clear that the models most 
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likely to avoid producing too much kSZ power will be the 
narrowest profiles consistent with the SN data: narrower 
profiles both produce less kSZ power when baryon sup- 
pression is ignored, and also result in greater baryon sup- 
pression, since the relevant length scales are the smallest. 



VI. IMPORTANCE OF THE CONSISTENT LTB 
FRAMEWORK 

At first glance, the kSZ power values in the 2j,-fij£ c 
plane presented in figure 6 appear remarkably similar 
to those of ZS10, considering the fully relativistic treat- 
ment of the background and the detailed treatment of 
the k scales performed here. However, note that for 
generic void profiles there is considerable arbitrariness 
in attempting to define an "edge" , with corresponding 
rcdshift zl- Therefore, we only expect plots of this kind 
for different profiles to agree up to a scaling of the Zl axis 
by a factor of roughly unity. Considering the strong sen- 
sitivity of kSZ power to Zl exhibited in figure 6, this arbi- 
trariness means that we can only say that our kSZ values 
are in agreement with those of ZS10 to within roughly 
an order of magnitude. This point is particularly rele- 
vant considering how different our smooth, polynomial 
fiducial LTB profile is from the "step function" Hubble 
bubble profile used in ZS10. 

While we cannot state how precisely our kSZ power 
values agree with those of ZS10, we can directly com- 
pare our kSZ values for models which satisfy the SN 
constraints. Superimposed on figure 6 are the 1, 2, and 
3cr confidence levels from the Union2 SN data calculated 
in section IV. Comparing these contours with those in 
ZS10, we can see that the kSZ power values for models 
which satisfy the SN constraints at 2a are roughly one 
tenth those found in ZS10. A related observation is that 
ZS10 emote a minimum \ 2 of 605.4 for their void models 
with respect to the 557 Union2 SNe, which represents a 
much worse fit than our best-fit value of x 2 = 539. We 
believe that the poorness of fit of the SZ10 profiles to the 
SNe accounts for most of the roughly tenfold difference 
in our results. This illustrates the importance of per- 
forming the kSZ calculations within a consistent, LTB 
framework, with smooth radial profiles which provide a 
good fit to the SNe. To obtain meaningful estimates of 
the kSZ power in void models, we cannot rely on void 
profiles which provide a poor fit to the SNe. 

We can reinforce this point by calculating a quantitity 
related to, but distinct from, the kSZ power, namely the 
Compton y-distortion of the CMB frequency spectrum. 
This distortion arises from the scattering of CMB pho- 
tons from inside our past light cone into our line of sight, 
and has been used to provide constraints on void mod- 
els [19, 58, 59]. Like the kSZ effect, the y-distortion also 
depends on the dipole /3(z) along our past light cone, but 
is a background-level effect, so is independent of the mat- 
ter power spectrum. In the single-scattering and linear 
approximations, and when the dipole anisotropy domi- 
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FIG. 9: Solid contours are the {/-distortion, as a function of 
central matter density parameter, f2^ c o, and a measure of the 
width of the void, zl- The contour values are log 10 y, calcu- 
lated using (40), and the heavy contour represents the COBE 
upper limit of y < 1.5 x 1CP 5 [60]. Results are presented for 
the fiducial LTB profile with ho = 0.71. Dotted contours are 
the 1, 2, and 3cr confidence levels from the SN data, which 
used the HST prior of h = 0.72 ± 0.08. 



nates, the y-distortion can be written as [19] 

y = To fj dr TzTr P{r)2 = To [J drF ^ r ^ ( 4 °) 

where r ro is the radial coordinate of (the assumed abrupt) 
reionization. 

The similarity of expression (40) to expression (30) for 
the kSZ power in the Limber approximation means that 
we can readily calculate the y-distortion for our fiducial 
profile. In figure 9 we plot the y-distortion in the zl-VL 1 ^ q 
plane for the fiducial LTB profile. The heavy contour 
indicates the 2a upper limit of y < 1.5 x 10~ 5 from the 
COBE satellite [60]. Superimposed on the figure are our 
1,2, and 3a confidence levels from the Union2 SN data. 
It is apparent that almost all of the models allowed by 
the SN data have y-distortion values below the COBE 
limit, and hence the y-distortion provides no significant 
constraint on void models. 

This result is in sharp contrast to [59] , who found that 
almost all models allowed by SN data were above the 
COBE limit, and hence ruled out. Importantly, the void 
models used in [59] were the same Hubble bubble mod- 
els as were studied in ZS10. The simplicity of the y- 
distortion calculation means that it is free of the ambigu- 
ities of harmonic decomposition, k scales, and small-scale 
baryonic astrophysics which plague the kSZ calculation. 
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Therefore, this result demonstrates again that the treat- 
ment of the background spacetime, and, in particular, 
obtaining good fits to the SN data, is crucial in obtaining 
meaningful constraints on void models for acceleration. 

VII. CONCLUSIONS 

Our results indicate that void models which satisfy 
local constraints on the Hubble rate and fit the SN 
observations predict considerably higher kSZ power for 
i ~ 2000-3000 than the total anisotropy power observed 
at those scales by SPT and ACT. Thus all local void mod- 
els would appear to be ruled out. To evade these high-£ 
constraints, the kSZ power would need to be reduced by 
a factor of roughly 30 for the narrowest voids allowed 
by the SN data, with our conservative interpretation of 
the power measurements. We have highlighted several 
caveats and ambiguities that render our calculations un- 
certain and could in principle lead to such a reduction. 
By way of summary, we list them here: 

(1) We have needed a harmonic decomposition of the 
matter fluctuations, although it is unclear how this 
should be done at late times on an LTB back- 
ground. 

(2) We have attempted to match the late time mat- 
ter power in the LTB model to that in standard 
AFLRW via the physically correct prescription, but 
we have found that it is not possible to do this 
without ambiguity in the k scales and hence in the 
matter power. 

(3) We have assumed that the linear kSZ contribution 
J Sf3dr vanishes in LTB models as it does in FLRW, 
although this would need to be checked explicitly 
with a proper LTB perturbation formalism. 

(4) The kSZ power is very sensitive to the parameters 
<78 and /b, but their local values are very uncertain, 
both theoretically and observationally. 

(5) We have shown that the kSZ power for large I is 
sourccd mainly by very small scales, between 1 and 
0.1 Mpc. The baryon power is expected to be sup- 
pressed significantly relative to total matter power 
on these scales, reducing the kSZ power by an un- 
certain amount. 

(6) The strong sensitivity of kSZ power to void width 
means that model dependence in the SN data, or 
simply changes to the void radial profile, may sub- 
stantially affect the kSZ power. 

(7) It is likely that the dipole values along the light 
cone, and hence the kSZ power, could be reduced 
by choosing an appropriate isocurvature profile be- 
tween radiation and matter at early times [25, 55], 
although it is likely that this would entail fine tun- 
ing. 



(8) It is also possible that including a significant late- 
time LTB decaying mode could evade the kSZ con- 
straints, although this would entail a drastic depar- 
ture from the standard view of the early Universe. 

To determine the importance of the first three points 
listed above would require a proper treatment of pertur- 
bations on LTB backgrounds. To reduce the uncertainty 
due to points (4) and (5) would require a direct measure- 
ment of the baryon power spectrum on ~l/3 Mpc scales, 
as well as a measurement of the baryon fraction, both at 
rcdshifts z ~ 0.5. It is clear, though, that a void must be 
as narrow as possible if it is to minimize the kSZ power. 

Our result that void models with the fiducial parame- 
ters predict, at 2a confidence, at least a factor of roughly 
30 larger kSZ power than that observed is considerably 
weaker than the corresponding result of ZS10, who found 
a discrepancy of more than three orders of magnitude for 
voids that fit the SNc at 2a confidence. The difference 
between our results is presumably mainly due to our rig- 
orous LTB treatment of the background together with 
our careful treatment of wave numbers. Our more con- 
servative interpretation of the SPT and ACT power mea- 
surements also contributes to the difference in our results. 
This conclusion is reinforced by the y-distortion values 
for our models, which are strongly at odds with those 
of [59], who also used a simple Hubble bubble model. 
Again, this emphasizes the importance of using smooth 
radial profiles which fit the SN data well in order to ob- 
tain meaningful constraints on void models. 

As we stressed in the Introduction, if a void model is 
to generate the observed primary CMB anisotropics, it 
must have a local Hubble rate low enough (ho ~ 0.45) to 
rule out all such models [16, 19, 20]. This conclusion is 
independent of the poorly understood details of the IC's 
or evolution of perturbations on LTB backgrounds, and 
it is insensitive to the uncertainties of small-scale bary- 
onic astrophysics as well as to the values of locally poorly 
determined parameters such as as and fb- For growing 
mode LTB profiles, the only chance to avoid this CMB 
constraint is possibly by substantially altering the pri- 
mordial spectrum from scale invariance [19]. But it is 
not enough, of course, to simply introduce a significant 
non-scale-invariant feature: that feature must be just 
so contrived that together with the local void profile it 
mimics the near-scale-invariant spectrum expected in the 
standard AFLRW model! In [44], e.g., a five-parameter 
modification of the primordial spectrum was required in 
order to match only the WMAP CMB data (no small- 
scale CMB data was used), with the somewhat improved 
local Hubble rate of ho ~ 0.60. The effect of an early 
radiation inhomogencity [55] should be negligible at late 
times [19], and hence it appears unlikely that this could 
provide a loophole. 

The superb sensitivity of the linear kSZ approach goes 
a very substantial way towards closing the rather con- 
trived loophole of modifying the primordial spectrum for 
void models. However, the various ambiguities in its ap- 
plication to void models should be considerably less im- 
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portant in a related but distinct application: constrain- 
ing departures from homogeneity within the standard 
AFLRW framework. That is, assuming that the back- 
ground is near-FLRW with the usual fraction of dark 
energy, and that the parameters erg and /b are as indi- 
cated by the CMB, we can use the linear kSZ effect to 
constrain small departures from homogeneity, as done in 
ZS10. Also, at greater redshifts, the length scales sourc- 
ing the kSZ at fixed I will be larger, reducing the un- 
certainty in baryonic astrophysics. Thus the linear kSZ 
effect should bring us closer to observationally confirming 
the long-held belief in cosmological homogeneity. 
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Appendix: The Lemaitre-Tolman-Bondi solution 

We model the Universe at background level as a spher- 
ically symmetric distribution of pressureless matter, with 
the observer at the centre. For this spacetime, Einstein's 
equations can be solved exactly, resulting in the LTB so- 
lution. The metric can be written 



ds 2 



-dt 2 



Y' 



1 - K 



dr 2 + Y 2 dn 2 



(A.l) 



where a prime denotes the derivative with respect to co- 
moving radial coordinate r, and t is the proper time along 
the comoving worldlincs. The function K = K(r) is ar- 
bitrary, and the areal radius Y = Y(t,r) is given para- 
metrically by 



f M 
~K 



(1 — cosh ?y) 



Y = < 



M 
~K 



-(1 — cos 77) 



(f) v V<.>- 



K < 0, 
< K < 1, 

K = 0, 



(A.2) 



M 



^§72 fa " smr ?) 



(sinluy — 77) K < 0, 



(A.3) 



< K < 1. 



This exact solution reveals the presence of a second free 
radial function, t& = ie (r) , which is known as the "bang 
time" function, since the cosmological singularity occurs 



at t = teM- Throughout this work we assume a homo- 
geneous big bang and set is = 0. In this case the LTB 
solution contains no decaying mode [17, 61], and hence 
we retain the standard inflationary picture of an essen- 
tially homogeneous early (but post-inflation) Universe. 
There is also a third free radial function, M(r), which 
we set to M(r) = r 3 without loss of generality as our 
gauge condition. 

The radial and transverse comoving expansion rates 
are given by Hh = Y'/Y' and H± = Y/Y, respec- 
tively, where the overdot denotes the derivative with re- 
spect to t. We define the local density parameter by 
fij^ c = 247rGp m /0 2 , where p m is the total matter density 
and 9 the comoving volume expansion. At the centre 
of symmetry in an arbitrary LTB spacetime, we have 
H\\ = H± = H. The angular diameter distance from the 
centre to rcdshift z is simply c?a = Y(z). 

The LTB spacetime is completely determined once the 
single free radial profile K(r) is specified. Based on our 
previous experience with a wide range of LTB profiles [16, 
19], we chose for this study a simple profile which depends 
on only two parameters, a width and a depth. Explicitly, 
we chose 



K{r) 



9 



KoL 
5 r 



r < L, 



r > L. 



(A.4) 

The profile width and depth are determined by the pa- 
rameters L and respectively, but it is normally more 
physically relevant to specify the parameter rij^ c = 

^m C ( z = 0) a t the observation point to characterize 
the depth, and to specify the redshift zl = z(r = L) 
to determine the width. Finally, the observation point 
is determined by the corresponding proper time coordi- 
nate, ten or i more observationally relevantly, the Hubble 
rate today at the centre, Hq = H(t<y). We will refer to 
the general profile, (A.4), with unspecified Ko, L, and 
Hq, as the fiducial profile. The model with the specific 
values ri^ c = 0.2, zl = 0.5, and ho = 0.71 (where 
Hq = 100 ho kms" 1 Mpc -1 ) will be referred to as the 
fiducial model; as shown in section IV, this model is a 
very good fit to the supernova data. In all of our calcu- 
lations we chose the values /b = 0.168 and as = 0.81 for 
the baryon fraction and matter amplitude today, as CMB 
observations suggest [32] in the context of FLRW mod- 
els. In our numerical code, we check for the presence of 
shell crossing singularities along the past light cone, and 
exclude any models which exhibit them, since the LTB 
solution is not valid in such Further details and a 

discussion of our numerical implementation are provided 
in [19]. 
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